Tunability of wire-grid metamaterial immersed into nematic liquid crystal 
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We propose electrically tunable hybrid metamaterial consisting of special wire grid immersed into 
nematic liquid crystal. The plasma-like permittivity of the structure can be substantially varied due 
to switching of the liquid crystal alignment by external voltages applied to the wires. Depending on 
the scale of the structure, the effect is available for both microwave and optical frequency ranges. 



Electromagnetic properties of periodic arrays of metal 
wires, stripes and finely structured particles have at- 
tracted increasing attention in recent years. The variety 
of emerging phenomena extends from the negative refrac- 
tion and subwavelength light constraining [1[ to extraor- 
dinary transmission [2[ , which all have been attributed to 
regular metallic structures with periods smaller or much 
smaller than the wavelength of electromagnetic radiation. 
As a rule, the unusual properties are characterized by a 
considerable or even anomalously strong frequency dis- 
persion. Undoubtedly, this opens prospective possibili- 
ties for tuning and switching, since a small modification 
of the system may result in a substantial change of its 
electromagnetic response at a given signal frequency. 

Liquid crystals (LC) inherently have all the necessary 
characteristics to provide tunability to subwavelength 
metal arrays. By switching the orientation of LC medium 
between metal elements or in metal cavities one can trig- 
ger substantial changes in the overall properties. Impor- 
tantly, LCs are transparent and significantly anisotropic 
in both microwave and optical frequency ranges, in which 
subwavelength metal arrays are usually operating |3|, [4[ . 
Finally, the metal constituents can be readily used as 
electrodes for applying switching voltages. 

In the past decade, various types of metamaterials have 
received significant attention due to negative values of 
permeability and permittivity (see e.g. the recent reviews 
H, Q and refs. therein). Certain methods are known 
how to provide tunability to metamaterials at microwave 
frequencies by inserting electronic semiconductor compo- 
nents 0? H, Q. However, it still remains challenging to 
achieve this at optical frequencies. In this Letter we de- 
scribe how nematic LC environment can provide tunabil- 
ity to the plasma frequency of wire-grid metamaterial. 

Wire lattices exhibit plasma-like permittivity for elec- 
tromagnetic radiation with wavelengths much larger than 
the lattice period UjJ. In particular, for electro- 
magnetic waves traveling normally to the grid of infinite 
parallel wires and polarized along them (xy-plane of in- 
cidence and z-polarization of electric field in Fig. [1]), the 
employed component of the effective local permittivity is 
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Here e^tz is the zz-component of the permittivity tensor 
of host medium surrounding the wires, and the parameter 
Q is determined by the grid geometry. 




FIG. 1: (Color online) Schematic of switching of nematic LC 
surrounding grid of wires. Note: the size of LC molecules is 
extremely exaggerated. 



The so-called plasma-frequency, at which the permit- 
tivity (pQ) changes its sign, is ujq = Q/^e^J) 1 / 2 . The grid 
on its own is transparent at frequencies higher than coo 
and reflecting at lower frequencies. Being accompanied 
by a 3D array of split-ring resonators providing negative 
permeability, the structure is a transparent left handed 
medium for oj < uoo and reflects waves with oj > uoq. 

The presence of ezz m Eq. (pQ) suggests an easy way 
of tuning the metamaterial by varying the host medium 
permittivity. Nematic LCs provide a natural opportunity 
to realize this. Anisotropic LC permittivity axes follow 
the orientation of the director n: 



(cj) = e±(u) Sij +s a (uj) mrij. 
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Typically, elongated molecules of nematic LC tend to 
align along interfaces due to surface anchoring. The re- 
sulting LC alignment along the wires is shown in Fig {]] on 
the left. If external voltages are applied to the neighbor- 
ing wires, this will load them with charges of different 
signs. The arising static electric field in the xy-plane 
will force the LC molecules to orient perpendicularly to 
the wires as shown in Fig IT] on the right. As a result, 
£zz switches from (e± + e a ) to e± providing the relative 
switching range of uoq 
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The dielectric anisotropy of LC is known to be at least of 
the order of several dozens of percent for both microwaves 
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FIG. 2: Grid capacitance per pair of wires in vacuum per 
length versus the ratio of wire diameter to lattice constant. 

and light 0, 0]- Therefore, the proposed design provides 
a possibility of shifting the plasma- frequency by 10-20%. 

Let us find out the voltage to be applied to the grid 
for efficient switching. Consider the quadratic lattice of 
cylindrical wires as shown in Fig. [H Assume that a volt- 
age V applied across the nearest neighbors loads the wires 
with charge per length densities ±Q. The charges are dis- 
tributed over the surfaces of thin wires. We neglect the 
angular inhomogeneity of this distribution and assume 
the electric field of a wire to be the same as that of a 
charged wire axis. We also presume that the static di- 
electric anisotropy of the nematic s s * — s a (0) is small and 
the electric field pattern is not perturbed by the inhomo- 
geneously oriented LC. 

A single wire stretched along the line x = 0, y = and 
loaded with the charge Q per unit length produces the 
electric field potential ipi(p) = — 27 r^ £st where 

p = \/ x 2 + y 2 , e 8 t = £_l(0) is the static permittivity of 
the LC, and the zero of the potential is assigned to the 
wire surface. In the infinite 2D lattice of alternatingly 
charged wires, the total potential outside of the wires is 

n,ra=oo 

p(r)= £ (-i) n+ "Vi Gw), (4) 

n,m= — oo 



where p nm = \/ (x — na) 2 + (y — ma) 2 and the coordi- 
nate origin lies on the axis of a positively charged wire. 

We have calculated (|4]) for finite latices. In spite of the 
absence of true mathematical convergence of the series, 
the potential distribution establishes already for a lattice 
of 10 x 10 wires. Further increase of the number of wires 
does not affect the potential in the middle provided that 
the overall electric neutrality is preserved. 

Applied voltage equals the difference of the potential 
on positively and negatively charged wires, V = (p+ — (f-. 
The capacitance per length of a grid of large even number 
N of wires is C = QN/(2V). It is convenient to intro- 
duce the capacitance per pair of wires in the vacuum 
per length, c\ = 2C/(Ne s t), which naturally appears in 
the relation Q = e st ciV. Apparently, c\ is a function 
of two geometric parameters, a and d, and depends on 
their ratio d/a. The calculated values of c\ are presented 
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FIG. 3: Lines of electric filed produced by 2D grid of alter- 
natingly charged wires. The wire diameter is set to be 0.1 of 
the lattice period 



in Fig. [21 Note that c\ tends to infinity at d ~ 1.08a 
instead of d = a, when the neighboring wire surfaces 
come in contact. This is the consequence of neglecting 
the anisotropy of the charge distribution over wire sur- 
face. We see that even for thick wires our approximation 
is reasonably good. 

Typical calculated electric field pattern is shown in 
Fig. 03 The field amplitude is maximum at the wire sur- 
faces (where the lines are condensed) and vanishes at the 
middle points of the diagonal lines connecting the nearest 
wires with the same sign of charge. 

Now we consider the effect of this field on the nematic. 
If we totally neglect the LC elasticity, the anisotropic 
molecules should align along the field lines. The director 
n then lies strictly in the xy-plane and the switching is 
perfect. The finite elasticity of the LC disturbs the ideal 
picture and the bending pattern of the director differs 
from the field lines pattern. However, this is of minor 
importance for the ei h J switching since the director still 
stays within the x?/-plane. 

The critical effect of LC elasticity occurs in the vicin- 
ity of the wires. To ensure the LC orientation along the 
wires when the voltage is switched off, the anchoring of 
LC molecules at wire surfaces must be strong. When 
the voltage is switched on, the anchoring still forces the 
molecules at the surfaces to point along the z-axis. Ac- 
cordingly, the director has to rotate by an angle of tt/2 
within the transient layer from the wire surface to the LC 
bulk. Efficient switching occurs when this layer is thin 
compared to the scale of the structure (lattice constant), 
i.e., practically important is the limit of strong electric 
field and thin transient layer. 

In this limit, we can neglect the contributions from 
other wires to the field near a wire surface. In the cylin- 
drical coordinates, the only present p-component of the 
electric field equals E p (p) = 27r ^ £ t P~ l while the nematic 
director has two components n p (p) and n z (p). At the 
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FIG. 4: Spatial profiles of cos 2 determining modulation of 
LC permittivity e zz - Corresponding values of v 2 are shown. 



wire surface n p (d/2) = 0. 

The free energy of the LC can be presented as the sum 
F = Fk + Fe of the elastic deformation energy and the 
dielectric energy [3], where the latter reads 



Jd 



F E = -7re £ s a / p dp (E p n p ) 



(5) 



Here the upper limit R is large enough to assure n z (R) — 
and can be extended to infinity. The simplest form 
of the LC elastic energy is given by the so-called one 
constant approximation: 



Fk = § / dV [ (v ' n)2 + (v x n)2 l 
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Introducing the director polar angle as n p = sind, 
n z = cos 9 we obtain the compact form of the total energy 



F = ttK 



P dp 



(e') 2 + (i-v 2 y- 
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where the parameter v 2 = V 2 c 2 e s ^ / {A-k 2 Ksq) character- 
izes relative contribution of the voltage-driven term. 

Minimizing the functional Q yields the differential 
equation for 0(p): 



1 

P 



y + (v 2 - 1) 



sin 20 
2p 2 







(8) 



with the boundary conditions 0(d/2) = and 0(oo) 
7r/2. For v 2 > 1 the exact solution reads: 
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The spatial variation of the zz-component of dielectric 
tensor ([2]), is controlled by the factor n 2 (p) = cos 2 6{p). 
The profiles of the latter are given in Fig.[Hfor several val- 
ues of v 2 larger than unity. It is seen that the voltage has 
to exceed considerably the critical value to provide the 
LC reorientation within thin transient layer. For v 2 > 2 
the layer thickness becomes comparable with the wire 
diameter. For estimates, we take a nematic LC elastic 
modulus K c± 10 _11 N, static permittivity e st ~ 10 and 
static dielectric anisotropy £ s * ~ 1. We also set the wire 
diameter to be ten times smaller than the lattice con- 
stant, which according to Fig. [2] yields c\ c± 10 _11 F/m. 
Then appropriate voltage is estimated as 



(10) 



Remarkably, this moderate voltage is independent of the 
scale of the grid. 

One notes that the studied LC switching geometry dif- 
fers qualitatively from the conventional one within fiat 
electrooptic cells. In our case, the LC bulk is directly 
aligned and reoriented by the immersed wires. We be- 
lieve this to be an additional advantage for microwave 
applications, where the controlling of millimeter (or even 
centimeter) thick LC samples by cell surfaces is practi- 
cally impossible. 

In summary, we have demonstrated that the plasma 
frequency of wire-grid metamaterial immersed into ne- 
matic LC can be efficiently tuned by 10-20% due to the 
switching of LC alignment. Our estimates show that the 
necessary voltage applied to the wires is of the order of 
several Volts. The main conclusions are valid for both 
microwave and optical frequency ranges. 
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